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Abstract 

We present a computer simulation model that is a one-to-one copy of an experimental realization of Wheeler's delayed choice 
experiment that employs a single photon source and a Mach-Zehnder interferometer composed of a 50/50 input beam splitter 
and a variable output beam splitter with adjustable reflection coeflicient R (V. Jacques et a]., Phys. Rev. Lett. 100, 220402 
(2008)). For < < 0.5, experimentally measured values of the interference visibility V and the path distinguishability 
D, a parameter quantifying the which-path information WPI, are found to fulflll the complementary relation + < 1, 
thereby allowing to obtain partial WPI while keeping interference with limited visibility. The simulation model that is solely 
based on experimental facts, that satisfles Einstein's criterion of local causality and that does not rely on any concept of 
quantum theory or of probability theory, reproduces quantitatively the averages calculated from quantum theory. Our results 
prove that it is possible to give a particle-only description of the experiment, that one can have full WPI even if D = 0, 
V = 1 and therefore that the relation + < 1 cannot be regarded as quantifying the notion of complementarity. 
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1. Introduction 

Particle- wave duality, a concept of quantum theory, 
attributes to photons the properties of both wave and 
particle behavior depending upon the circumstances of 
the experiment [1]. The particle behavior of photons 
has been shown in an experiment composed of a single 
beam splitter (BS) and a source emitting single pho- 
tons and pairs of photons [2]. The wave character has 
been demonstrated in a single-photon Mach-Zehnder 
interferometer (MZI) experiment [2]. In 1978, Wheeler 
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proposed a gedanken experiment [3], a variation on 
Young's double slit experiment, in which the decision 
to observe wave or particle behavior is made after the 
photon has passed the slits. The pictorial description 
of this experiment defies common sense: The behavior 
of the photon in the past is said to be changing from a 
particle to a wave or vice versa. 

Recently, Jacques et al. reported on an experimen- 
tal realization of Wheeler's delayed choice experiment 
using a single photon source and a MZI composed of 
a 50/50 input BS and a variable output BS with ad- 
justable reflection coefficient R [4], a modification of 
the experiment presented in [5] in which two 50/50 BSs 
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Fig. 1. Schematic diagram of the experimental setup for Wheeler's delayed-choice gedanken experiment [4]. PBS: Polarizing 
beam splitter; HWP: Half-wave plate; EOM: electro-optic modulator; RNG: Random number generator; WP: Wollaston 
prism; P,S: Polarization state of the photon; ^: Phase shift between paths and 1; Dq, D\: Detectors. 



were used. A schematic picture of the experimental set- 
up is shown in Fig. 1. The reflection coefficient R of the 
variable beam splitter (BSontpnt) can be controlled by a 
voltage applied to an electro-optic modulator (EOM), 
making it act as a variable wave plate. This can be done 
after each photon has entered the MZI. The phase-shift 
$ between the two arms of the MZI is varied by tilt- 
ing the polarizing beam splitter (PBS) of the variable 
output BS. For < < 0.5 measured values of the in- 
terference visibility V [6] and the path distinguishabil- 
ity D [4], a parameter that quantifies the which-path 
information (WPI), were found to fulfill the comple- 
mentary relation -\- < 1. The extreme situations 
(y = 0, D = 1) and {V = 1, D = 0), obtained for 
R = and R = 0.5, give full and no WPI, associated 
with particlelike and wavelike behavior, respectively. 
By choosing < R < 0.5 Jacques et al. claim to have 
obtained partial WPI while keeping interference with 
limited visibility [4], thereby having accomplished an 
affirmative delayed choice test of complementarity or 
wave-particle duality as it is often phrased. 

Although the detection events (detector "clicks") are 
the only experimental facts, the pictorial description of 
Jacques et al. [4,5] is as follows: Linearly polarized sin- 
gle photons are sent through a 50/50 PBS (BS^nput), 



spatially separating photons with S polarization (path 
0) and P polarization (path 1) with equal frequencies. 
After the photon has passed BSinpnt, but before the 
photon enters the variable output the decision to ap- 
ply a voltage to the EOM is made. The PBS of BS output 
merges the paths of the orthogonally polarized photons 
travelling paths and 1 of the MZI, but afterwards the 
photons can still be unambiguously identified by their 
polarizations. If no voltage is applied to the EOM then 
R — and the EOM can be regarded as doing nothing 
to the photons. Because the polarization eigenstates 
of the Wollaston prism correspond to the P and S po- 
larization of the photons travelling path an 1 of the 
MZI, each detection event registered by one of the two 
detectors Do or Di is associated with a specific path 
(path or 1, respectively). Both detectors register an 
equal amount of detection events, independent of the 
phase shift $ in the MZI. This experimental setting, 
corresponding to the open configuration of the MZI, 
clearly gives full WPI about the photon within the in- 
terferometer (particle behavior), characterized by D = 
1. In this case no interference effects are observed, cor- 
responding with a zero interference visibility (y = 0). 
When a voltage is applied to the EOM, then R ^ 
(see Eq. (2) in [4]) and the EOM acts as a wave plate 
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rotating the polarization of the incoming photon by an 
angle depending on R. The Wollaston prism partially 
recombines the polarization of the photons that have 
travelled along different optical paths with phase dif- 
ference $ (closed configuration), and interference ap- 
pears (y / 0), a result expected for a wave. The WPI 
is partially washed out, up to be totally erased when 

R = o.5 (D^o, y = 1). 

The outcome of delayed-choice experiments [4,5,7,8,9,10, 
that is the average results of many detection events, is 
in agreement with wave theory (Maxwell or quantum 
theory). However, the pictorial description explain- 
ing the experimental facts [5] defies common sense: 
The decision to apply a voltage to the EOM after 
the photon left BSinput but before it passes BSoutput, 
influences the behavior of the photon in the past and 
changes the representation of the photon from a par- 
ticle to a wave [5]. Although on one hand quantum 
theory can be used to describe the final outcome of 
this type of experiments (the average results of many 
detection events), on the other hand it does not de- 
scribe single events [1]. Therefore, it should not be a 
surprise that the application of concepts of quantum 
theory to the domain of individual events may lead 
to conclusions that are at odds with common sense. 
Although not applying this reasoning to describe this 
type of experiments could prevent us from making 
nonsensical conclusions, this unfortunately would not 
give us a single clue as how to explain the fact that in- 
dividual events are observed experimentally and, when 
collected over a sufficiently long time, yield averages 
that agree with quantum theory. 

Since no theory seems to exist that can give a sensi- 
cal description of the "whole" experiment, we adopted 
the idea to search for algorithms that could mimic 
(simulate) the detection events and experimental pro- 
cesses, including for example the random switching of 
the EOM for each photon sent into the interferometer. 
We moreover require that the algorithms used to sim- 
ulate the action of an optical element, such as a BS or 
wave plate, on a photon should be independent of the 
experimental setup. In other words, the algorithms to 
simulate an optical component should be the same for 
all identical optical components within the same exper- 
iment but also within a different experiment. Hence, 
first solving the Schrodinger equation for a given ex- 
perimental configuration and then simply generating 
events according to the resulting probability distribu- 



tion is not what we have in mind when we perform 
an event- by- event simulation of the experiment. Sim- 
ilarly, first calculating the quantum potential (which 
requires the solution of the Schrodinger equation) and 
then solving for the Bohm trajectories is a different 
kind of event- by- event simulation than the one we de- 
scribe in this paper. In this paper, the event-by-event 
simulation algorithm reproduces the results of quan- 
llj,^m theory, without first solving a wave equation. 

In this paper, we describe a model that, when im- 
plemented as a computer program, performs an event- 
by-event simulation of Wheeler's delayed-choice exper- 
iment. Every essential component of the laboratory ex- 
periment (PBS, EOM, HWP, Wollaston prism, detec- 
tor) has a counterpart in the algorithm. The data is an- 
alyzed by counting detection events, just like in the ex- 
periment [4,5]. The simulation model is solely based on 
experimental facts, satisfies Einstein's criterion of local 
causality and does not rely on any concept of quantum 
theory or of probability theory. Nevertheless, our simu- 
lation model reproduces the averages obtained from the 
quantum theoretical description of Wheeler's delayed 
choice experiment but as our approach does not rely 
on concepts of quantum theory and gives a description 
on the level of individual events, it provides a descrip- 
tion of the experimental facts that does not defy com- 
mon sense. In a pictorial description of our simulation 
model, we may speak about "photons" generating the 
detection events. However, these so-called photons, as 
we will call them in the sequel, are elements of a model 
or theory for the real laboratory experiment only. The 
experimental facts are the settings of the various ap- 
paratuses and the detection events. What happens in 
between activating the source and the registration of 
the detection events is not measured and is therefore 
not known. Although we always have full WPI of the 
photons in the closed configuration of the interferom- 
eter (we can always track the photons during the sim- 
ulation) , the photons build up an interference pattern 
at the detector. The appearance of an interference pat- 
tern is commonly considered to be characteristic for a 
wave. In this paper, we demonstrate that, as in exper- 
iment, it can also be build up by many photons. These 
photons have full WPI, never directly communicate 
with each other and arrive one by one at a detector. 

The work described in this paper elaborates on the 
work described in [12] to simulate the experiment re- 
ported in [5]. The simulation model is built on ear- 
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lier work [13,14,15,16,17,18,19,20,21,22] that demon- 
strated that it may be possible to simulate quantum 
phenomena on the level of individual events without 
invoking a single concept of quantum theory. Specif- 
ically, we have demonstrated that locally-connected 
networks of processing units with a primitive learn- 
ing capability can simulate event-by-event, the single- 
photon beam splitter and MZI experiments of Grang- 
ier et al. [2] and Einstein-Podolsky- Rosen experiments 
with photons [23,24,25]. Furthermore, we have shown 
that this approach can be generalized to simulate uni- 
versal quantum computation by an event- by- event pro- 
cess [15,16]. Our event-by-event simulation approach 
rigorously satisfies Einstein's criterion of local causal- 
ity and builds up the final outcome that agrees with 
quantum theory event-by-event, as observed in real ex- 
periments. 



2. Simulation model 

The simulation algorithm can be viewed as a 
message-processing and message-passing process: It 
routes messengers one by one through a network of 
units that process messages. The messengers may be 
regarded as "particles" . These messengers carry a mes- 
sage which contains information about for example 
the relative time the particle traveled, its polarization, 
its color, its velocity, and so on. In other words, the 
message represents a so-called variable property of the 
particle that can be manipulated and measured given 
particular experimental setttings. The components 
of the experimental setup such as the BS, the wave 
plates, the Wollaston prism and so on are so-called 
processing units that interpret and manipulate the 
messages carried by the particles. These processing 
units are put in a network that represents the complete 
experimental setup. Since at any given time there is 
only one messenger being routed through the whole 
network, there is no direct communication between 
the messengers. The only form of communication is 
through the processing units when the messengers are 
routed through the network. The model satisfies the 
intuitive notion of local causality. 

In general, processing units consist of an input stage, 
a transformation stage, an output stage and have an 
internal vector representing their internal state. The 



input (output) stage may have several channels at 
(through) which messengers arrive (leave). Some pro- 
cessing units are simpler in the sense that the input 
stage is not necessary for the proper functioning of the 
device. As a messenger arrives at an input channel of 
a processing unit, the input stage updates its inter- 
nal vector, and sends the message, represented by a 
vector, together with its internal vector to the trans- 
formation stage that implements the operation of the 
particular device. Then, a new message is sent to the 
output stage, using a pseudo-random number to se- 
lect the output channel through which the messenger 
will leave the unit. We use pseudo-random numbers 
to mimic the apparent unpredictability of the experi- 
mental data only. The use of pseudo-random numbers 
is merely convenient, not essential. 

In the experimental realization of Wheeler's delayed 
choice experiment by Jacques et al. [4] linearly polar- 
ized single photons are sent through a PBS that to- 
gether with a second, movable, variable output PBS 
with adjustable reflectivity R forms an interferometer 
(see Fig. 1). The basic idea now is that we have to con- 
struct a model for the messengers representing the pho- 
tons and for the processing units representing the opti- 
cal components in the experimental setup. We require 
that the processing units for identical optical compo- 
nents should be reusable within the same and whitin 
different experiments. The network of processing units 
is a one-to-one image of the experimental setup [4] and 
is shown in Fig. 1. In what follows we describe some 
elements of our model in more detail. Additional infor- 
mation can be found in [12]. 

2.1. Messengers 

In a pictorial description of the experiment the pho- 
tons can be regarded as particles playing the role of 
messengers. Each messenger carries a (variable) mes- 
sage which contains information about its phase and 
polarization. The phase combines information about 
the frequency of the light source and the time that par- 
ticles need to travel a given path. However, no explicit 
information about distances and frequencies is required 
since we can always work with relative phases. 

The information carried by the messenger can be 
represented by a six-dimensional unit vector y^^^ = 
(cos ipk,n, sin ip^^^, cos V'fc^^, sin cos ^k,n, sin ^k,n)- 
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Fig. 2. Diagram of a DLM-based processing unit that performs an event-based simulation of a polarizing beam splitter 
(PBS). The solid lines represent the input and output channels of the PBS. The presence of a message is indicated by an 
arrow on the corresponding channel line. The dashed lines indicate the data flow within the PBS. 



The superscript H (V) refers to the horizontal (verti- 
cal) component of the polarization and ipk,n^ i^k,n^ 
^k,n represent the phases and polarization of the pho- 
ton, respectively. It is evident that the representation 
used here maps one-to-one to the plane- wave descrip- 
tion of a classical electromagnetic field [6] , except that 
we assign these properties to each individual photon, 
not to a wave. The subscript n > numbers the con- 
secutive messages and k = 0,1 labels the channel of 
the PBS at which the message arrives (see below). 

Since in this paper we will demonstrate explicitly 
that in our model photons always have full WPI even 
if interference is observed, we give the messengers one 
extra label, the path label having the value or 1. The 
information contained in this label is not accessible in 
the experiment [4] . We only use it to track the photons 
in the network. The path label is set in the input BS 
and remains unchanged until detection. Therefore we 
do not consider this label in the description of the pro- 
cessing units but take it into account when we detect 
the photons. 



2.2. Polarizing beam splitter 

From classical electrodynamics we know that if an 
electric field is applied to a dielectric material the ma- 
terial becomes polarized [6]. The polarization P(k, t) 
is given by 

P(k,t)= f xiKu)E{k,t-u)du, (1) 
Jo 

where E(k, t) denotes the electric field vector, k is the 
wave vector, and x is the linear response function [6]. 
From Eq. (1) it is evident that the dielectric material 
shows some kind of memory effect because the response 
(the polarization) of the material to the applied electric 
field is a function of both present and past values of 
the electric field. We use this kind of memory effect in 
our algorithm to model the PBS. 

The processor that performs the event- by- event sim- 
ulation of a PBS is depicted in Fig. 2. It consists of an 
input stage, a simple deterministic learning machine 
(DLM) [13,14,15,16], a transformation stage (T), an 
output stage (O) and has two input and two output 
channels labeled with k = 0,1. We now define the op- 
eration of each stage explicitly. 

- Input stage: The DLM receives a message on ei- 
ther input channel or 1, never on both chan- 
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nels simultaneously. The arrival of a message on 
channel (1) is named a (1) event. The input 
events are represented by the vectors Vn = (I5O) 
or Vn = (O5 1) if the nth event occurred on chan- 
nel or 1, respectively. The DLM has six internal 
registers Y|^^ — (C^^^^ ^k,n)^ ^k,n — (Ck,n^ ^k,n)^ 
Yj^n = {^k,n^ i^k.n) ^^nd ouc internal vector = 
), where xo,n + xi^n — 1 and > 

for A; = 0, 1 and all n. These seven two-dimensional 
vectors are labeled by the message number n be- 
cause their contents are updated every time the 
DLM receives a message. Note that the DLM stores 
information about the last message only. The infor- 
mation carried by earlier messages is overwritten by 
updating the internal registers. 

Upon receiving the (n + l)th input event, the 
DLM performs the following steps: It stores the 
first two elements of message yfc,n+i in its internal 
register Y^^^+i = {Ck^^^i, S^^^^i), the middle two 
elements of yfc,n+i in Y^^+i = {Ck,n+i, Sk,n+i)^ 
and the last two elements of yk n+i in n+i = 
n+i5 *^jr,n+i)- Then, it updates its internal vec- 
tor according to the rule [13] 

Xi,n+1 = OLXi^n + (1 " OL)^i,k', (2) 

where < a < 1 is a parameter that controls the 
learning process [13]. Note that by construction 

Xo,n^\ + a:i,n+l = 1, a:0,n+l > and Xi,n+1 > 0. 

From the solution of Eq. (2), 

n-l 

Xn = a'^xo + (1 - a) ^ OL^~'^~^Wj, (3) 

the correspondence to the expression for the polar- 
ization in classical electrodynamics Eq. (1) can be 
seen. The vector v plays the role of the electric field 
vector E and the internal vector x plays the role of 
the polarization P. Hence, one could say that the in- 
ternal vector X is the response of the PBS to the in- 
coming messages (photons) represented by the vec- 
tors V. Therefore the PBS "learns" so to speak from 
the information carried by the photons. The charac- 
teristics of the learning process depend on the pa- 
rameter OL (corresponding to the response function). 
Equation (2) is the simplest learning rule we could 
think of. If experimental measurements for a single 
PBS would require another maybe more complicated 
rule to simulate the experimental outcome then we 



could modify the learning rule but given the infor- 
mation we have right now Eq. (2) suffices. 
- Transformation stage: The second stage (T) of 
the DLM-based processor takes as input the data 
stored in the six internal registers ^k.n^i — 

(Ck,n-\-l^ ^k,n-\-l)^ = (Ck^,n-\-l^ ^X,n-\-l) ^ 

Yjr,n+i = (CCn+i'*5'j^n+i) ^ud in the internal 
vector Xn+i = (a:o,n+i, a:i,n+i) and combines the 
data into an eight-dimensional vector (see Fig. 2). 
Rewriting this vector as 
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(4) 

shows that the operation performed by the trans- 
formation stage T corresponds to the matrix-vector 
multiplication in the quantum theoretical descrip- 
tion of a PBS, namely 
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where (a^ , aX , , aY) and (bo ,bX ,bY) denote 
the input and output amplitudes of the photons with 
polarization H and V in the and 1 channels of a 
PBS, respectively. Note that in our simulation model 
there is no need to introduce the concept of a vac- 
uum field, a requirement in the quantum optical de- 
scription of a PBS. 
- Output stage: The final stage (O) sends the message 
w = {wo, wi,W2,W3, W4, W5) , where 
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^0 = C'^n+lC'^n+lV^0,n+l/^^^4, 
Wi = So^ri-\-lCo^n+ly/Xo,n+l/uW4, 
W2 = -SY^ri+lS[^n+lyyXl,n+l/uW5, 
'WS = CY,n+lS[^n+lVxi~^/uW5, 

(6) 

through output channel if > r where < r < 
1 is a uniform pseudo-random number. Otherwise, 
if It < r, the output stage sends through output 
channel 1 the message z = (2^0,^1,2:2,2:3,^4,^5) , 
where 

Z2 — — So^ri+lSo^n+l^Xo^n+l /VZ5, 
^3 — Co^n+lSo,n+lVxo~^/vZ5, 
Z4 = ^Zl+Zl/V, 

2:5 = \lzl+zl/v, 
v=^zl+zl+zl+zl. 

(7) 

Any other algorithm that selects the output channel in 
a systematic manner might be employed as well. This 
will change the order in which messages are being pro- 
cessed but the content of the messages will be left intact 
and the resulting averages do not change significantly. 

2.3. Remaining optical components 

The Wollaston prism is a PBS with one input channel 
and two output channels and is simulated as the PBS 
described earlier. 

In contrast to the PBS, the HWP and the EOM are 
passive devices. As can be seen from the wave mechan- 
ical description, a HWP does not only change the po- 
larization of the photon but also its phase [6] . 

When a voltage is applied to the EOM, R ^ (see 
Eq. (2) in [4]) and the EOM acts as a wave plate that 
rotates the polarization of the incoming photons by 



an angle depending on R. In the simulation a pseudo- 
random number is used to decide to apply a voltage to 
the EOM or not. Also here we use a pseudo-random 
number to mimic the experimental procedure to con- 
trol the EOM [4,5]. Any other (systematic) sequence 
to control the EOM can be used as well. 

2.4. Detection and data analysis procedure 

Detector Dq (Di) registers the output events at 
channel (1). During a run of N events, the algorithm 
generates the data set 

r{R) = {Xn,yn, An\n = 1, TV; $ = $i - , (8) 

where Xn = 0, 1 indicates which detector fired (Do or 
Di), yn = 0, 1 indicates through which arm of the MZI 
the messenger (photon) came that generated the detec- 
tion event (note that yn is only measured in the simula- 
tion, not in the experiment), and An = 0, 1 is a pseudo- 
random number that is chosen after the nth messenger 
has passed the first PBS, determining whether or not 
a voltage is applied to the EOM (hence whether the 
MZI configuration is open or closed). Note that in one 
run of N events a choice is made between no voltage 
or a particular voltage corresponding to a certain re- 
flectivity R of the output BS (see Eq. (2) in [4]). The 
angle $ denotes the phase shift between the two inter- 
ferometer arms. This phase shift is varied by applying 
a plane rotation on the phase of the particles entering 
channel of the second PBS. This corresponds to tilt- 
ing the second PBS in the laboratory experiment [4]. 
For each $ and MZI configuration the number of (1) 
output events A^o (Ni) is calculated. 



3. Simulation results 

The algorithm described above directly translates 
into a simple computer program that simulates the 
messenger routing in a network that contains all the 
optical components of the laboratory experiment [4]. 
Before the simulation starts we set xo = (xo,o,a:i,o) = 
(r, 1 — r) , where r is a uniform pseudo-random num- 
ber. In a similar way we use pseudo-random numbers 
to initialize Y^o, 'yX,o, Y^o, Yf^o, 'yY,o and Y^o- In 
this simulation, we send messengers to one input chan- 
nel of the input PBS only (see Fig. 1). The HWP in 
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Fig. 4. (Color online) Event-by-event simulation results of the interference visibility V for different values of R ((a)-(c)) and 
of V^, D'^ and V"^ + D'^ as a function of the EOM voltage (d). (a)-(c) Circles give the results for the normalized intensities 
Nq/N and Ni/N as a function of the phase shift A^o (Ni) denoting the number of events registered at detector Dq (Di), 
for (a) R = 0.43 (V ^ 0.98), (b) R = 0.05 (V ^ 0.45) and (c) R = (V = 0). For each value of the number of input 
events N = 10000. The number of detection events per data point is approximately the same as in experiment. Dashed lines 
represent the results of quantum theory (Malus law), (d) Circles give the simulation results. Lines represent the theoretical 
expectations obtained from Eqs. (2), (3) and (7) in [4] with (3 = 24^ and K = 2171^. 



^Soutput changes the phases and also interchanges the 
roles of channels and 1. Disregarding a few excep- 
tional events, the PBS in BSoutput generates messages 
in one of the channels only. For a fixed set of input 
parameters, each simulation takes a few seconds on a 
present-day PC. In all simulations, a = 0.99 [13]. 

We first demonstrate that our model yields full WPI 
of the photons. Figure 3 shows the number of detection 
events at Do as a function of $ for R = 0.5. The events 
generated by photons following path and path 1 of 
the MZI are counted separately. It is clear that the 
number of photons that followed path and path 1 is 
equal and that the total intensity in output channel 
obeys Malus law. Hence, although the photons have full 
WPI for all $ they can build an interference pattern 
by arriving one by one at a detector. 

Next, we calculate for = 0, 0.05, 0.43 [4] and for 
each phase shift $ and configuration (open or closed) 
of the MZI the number of events registered by the two 



detectors behind the output BS, just like in the experi- 
ment. Figure 4(a)- (c) depicts the interference visibility 
V. The simulation data quantitatively agree with the 
averages calculated from quantum theory and qualita- 
tively agree with experiment (see Fig. 3 in [4]). Calcu- 
lation of D as described in [4] gives the results for 
and V"^ shown in Fig. 4(d). Comparison with Fig. 4 in 
[4] shows excellent qualitative agreement. 



4. Conclusion 

In this paper, we have presented a simulation model 
that is solely based on experimental facts, that sat- 
isfies Einstein's criterion of local causality, that does 
not rely on any concept of quantum theory or of prob- 
ability theory, and that provides a description of the 
experimental observations in [4] on the level of indi- 
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Phase Shift 

Fig. 3. (Color online) Event-by-event simulation results of 
the interference visibility V for R = 0.5. Markers give the 
results for the normalized intensity Nq/N as a function 
of the phase shift <l>, A^o denoting the number of events 
registered at detector Dq. Circles (triangles) represent the 
detection events generated by photons that followed path 
(1) and squares represent the total number of detection 
events. For each value of the number of input events 
= 10000. The total number of detection events per data 
point (squares) is approximately the same as in experiment. 
The solid line represents the results of quantum theory. 

vidual events. In a pictorial description of our simula- 
tion model, we may speak about "photons" generating 
the detection events. In the simulation we can always 
track the photons, even in the closed configuration of 
the MZL The photons always have full WPI, never di- 
rectly communicate with each other, arrive one by one 
at a detector but nevertheless build up an interference 
pattern at the detector in the case of the closed config- 
uration of the MZL Hence, although for < i? < 0.5 
we find that < D <1 and + V'^ <1 with values 
for D and V in qualitative agreement with the experi- 
mental results, we always have access to full WPI, even 
in the case D = 0, V = 1. Our model thus provides a 
counter example for the fact that full WPI would cor- 
respond to D = 1. A further consequence is that the 
relation -\- < 1 cannot be regarded as quantify- 
ing the notion of complementarity: Our model allows a 
particle-only description for both the open and closed 
configuration of the MZL 
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